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Two-flavour Schwinger model with dynamical fermions in the Liischer 
formalism 

S. Elser * and B. Bunk t 

Institut fiir Physik, Humboldt-Universitat, Invalidenstr.llO, 10115 Berlin, Germany 

We report preliminary results for 2D massive QED with two flavours of Wilson fermions, using the Hermitean 
variant of Liischer's bosonization technique. The chiral condensate and meson masses are obtained. The simplicity 
of the model allows for high statistics simulations close to the chiral and continuum limit, both in the quenched 
approximation and with dynamical fermions. 



1. Schwinger model 

Studies of algorithms for dynamical fermions 
are very time-consuming, therefore we propose 
the massive Schwinger model (QED in two dimen- 
sions) as a low-cost laboratory. We start from the 
Euclidean continuum Lagrangean 
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with two degenerate flavours of fermions. For 
m = it is exactly soluble [|l| and has a known 
expansion in the mass|Q. Close to the chiral 
limit, we expect three light states ('pions'), mas- 
sive mesons and a non-zero fermion condensate. 

The lattice version is defined using compact 
link variables U^{x) G U{1) with the standard 
plaquette action Sg [U] and Wilson fermions with 
operator M. Integration over the fermion vari- 
ables results in a positive effective action for the 
gauge fields: 



Pcs[U] 



detM-^ e 



2. Liischer's local bosonic theory 

As an alternative to the Hybrid Monte Carlo 
algorithm, M. Liischer proposed a local bosonic 
formulation 1^. Let Q = cj^M = be scaled 
so that its eigenvalues are in [—1, 1] and P„(s) a 
polynomial of even degree n which approximates 
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1/s in (0,1]. Its roots Zk {k = l...n) come in 
complex conjugate pairs and determine = 
fJ-k + ivk {vk > 0). Then 

Fcff[C/] cx detg^e-^sl^l 

= det[g2p„(Q2)] [dg^ p„(g2)]-i ^-s,[u] 



with C(g2) EE det[Q^P„{Q^)] « 1 and n complex 
bosonic Dirac fields (pk- 

3. Implementation 

We chose as approximation polynomials P„(s) 
the Chebychev polynomials proposed by Bunk et 
al.Q. The convergence of P„(s) 1/s as n — > oo 
is exponential and uniform for s G [£,1]; introduc- 
ing a parameter e into the algorithm. 

The updating process consists of exact heat 
bath sweeps for the (/)'s and C/'sQ, followed by 
a number of over-relaxation iterations. We con- 
firmed, in our preliminary runs, that the reflec- 
tion sweeps for the ^'s and C/'s have to be com- 
bined in pairs, as was observed before Q. Finally, 
the approximation C « 1 has to be controlled. 
To this end we compute the lowest 8 eigenval- 
ues!^ Q^' them to estimate the change in 
C and apply a global Metropolis correction step 
with acceptance probability min [1,^]. This will 
be improved in the future. 

In order to avoid trouble with bad pseudo- 
random numbers we use Liischer's high-quality 
random number generatorlq]. 
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4. Observables 

We measured averages of local fermion bilin- 
ears with a noisy estimator scheme, using ran- 
dom spinors r]{x) = ±1 as sources for the CG 
solver. (V'TsV') {ipjfj,ip) were checked to van- 
ish within errors. The fermion condensate (ip^p) 
remains nonzero even in the chiral limit because 
Wilson fermions break chiral invariance explicitly. 
Our MC results are presented below. 

For the determination of meson masses opera- 
tors with various quantum numbers are defined: 

• flavour triplet - (ip^sT'ip) 'tt', {tpTip) 'ao', 

• flavour singlet - (^'75^/') (V'V') '/o' • 

Moreover, insertion of a 70, e.g. (jpjoj5T^p) for 
the TT, leads to alternative operators with the 
same quantum numbers in the rest frame. This 
exhausts the Dirac algebra in two dimensions. 

The calculation of temporal correlators r(At) 
involved point sources at randomly chosen posi- 
tions {x,t) and summation over spatial displace- 
ments y in the other time slice {y, t + At) to 
project out the zero momentum states. As to the 
flavour-singlet channels, the disconnected piece 
was subtracted with the aid of the noisy inver- 
sions performed earlier for the measurement of 
the condensate. 

5. Calculations 

5.1. Dynamical condensates 

The condensate with dynamical fermions was 
calculated on 16x16 lattices for (3 — 2.0, k = 
0.20 . . . 0.275 and on 16x32 lattices for /3 = 10.0, 
K = 0.20... 0.25. 

The Liischer algorithm required a small num- 
ber (n = 20 . . . 40) of boson fields only and al- 
lowed for very high values of e = 0.01. ..0.1. 
With moderate statistics (200 . . . 400 measure- 
ments), we reproduced the results obtained with 
a Hybrid Monte Carlo by the Graz group see 
Fig. 0. 

5.2. Quenched masses 

Meson masses in the quenched approxima- 
tion were obtained on 16x32 lattices for /3 — 
2.0, 6.0, 10.0, 20.0 and k = 0.20 . . . 0.275 perform- 
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Figure 1. Dynamical fermion condensate on a 
16x16 lattice at /3 = 2.0 



ing high-statistics runs with about 2000 indepen- 
dent measurements each. In case oi (3 — 6.0, Fig. 
1^ shows the clear signal for the pion mass decreas- 
ing as K, is increased, its alternative operator gives 
masses which agree within errors. For the 77 and 
its variant, consistency is also verified with larger 
errors. The higher states give noisy results. 

5.3. Dynamical masses 

With dynamical fermions, Irving et al.|l^] gave 
a mass of 0.369(3) for a 32x32 lattice with /3 = 
2.29, K — 0.26. We checked our implementation 
of the dynamical fermion update with n — 40 and 
obtained m — 0.377(4) from about 2000 measure- 
ments. 

More masses were determined on 16x32 lattices 
for (3 = 10.0, K = 0.20... 0.25 and are shown 
in Fig. |3[ We observe the same qualitative be- 
haviour as in the quenched case. 

6. Conclusions and outlook 

The two-flavour Schwinger model has physi- 
cal properties similar to QCD in four dimensions. 
It is much easier to simulate even with dynami- 
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Figure 2. Quenched meson masses on a 16x32 
lattice at /3 = 6.0 
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Figure 3. Dynamical meson masses on a 16x32 
lattice at /? = 10.0 



cal fermions and allows to determine observables 
with high precision. The preliminary results pre- 
sented above are encouraging. 

The next steps will be to include standard 
versions of preconditioning, a test of the non- 
hcrmitcan variant of the polynomial approxi- 
mation and a noisy estimator which makes the 
Metropolis acceptance step exact. 

As to the correlation functions, different noisy 
estimator schemes will be tested and their effi- 
ciency compared. 

Finally, the Schwinger model will be used to 
tune the boson algorithm and to determine how 
the computational cost scales in the continuum 
limit am^ with m-j^/mri and Lm^ kept fixed. 
Although larger lattices than those used so far 
will be necessary, the cost will be lower by orders 
of magnitude as compared to the case of QCD 
in four dimensions. This makes the Schwinger 
model a reasonable testing ground for dynamical 
fermion algorithms. 
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